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1 Introduction 



Super Yang-Mills theories (SYM) have attracted much interests as leading candidates for 
physics beyond the standard model. To understand their nonperturbative aspects, one might 
expect lattice simulation to be useful as it happens in the usual non-supersymmetric gauge theories. 
However, it is generally difficult to realize SYM on lattices and consequently detailed numerical 
simulations have been performed only in lower (less than four) dimensions (see pQ for recent 
simulations of 4d SYM, and [2] for a review of lattice supersymmetry) . 

SYM's are also important because it is expected to provide a nonperturbative formulation of 
superstring/M theory at large- iV [31 HJ [5j El [7] . In this context the relevant theories are SYM's 
in lower dimensional spacetime. In particular, (0 + l)-dimensional theory can be analyzed on 
computer by using the non-lattice technique [8] (in this case lattice simulations are also possible 
[9]) and a part of these conjectures has been confirmed by Monte-Carlo simulations in the strong 
coupling regime [101 fTT| 112] . Chern-Simons gauge theories are also relevant in the context of 
superstring/M theory. In fact recently supersymmetric Chern-Simons gauge theories in three 
dimensions have been proposed as a description of the theory of multiple M2-branes [13\ I14j . 
Since these theories describe membranes in their strong coupling regimes, numerical simulations 
appear to be a very useful tool. Unfortunately it turns out to be a very difficult task to study 
either Chern-Simons or supersymmetry on a lattice. 

At large- iV, it is possible to circumvent the lattice-SUSY problem by using the Eguchi-Kawai 
equivalence [15]. This construction guarantees that the large- iV gauge theories are equivalent to 
the lower dimensional matrix models if a certain condition is satisfied. Recently this equivalence 
has been used to formulate 4d M = 4 SYM compactified on S 3 [16] In this construction, the 
BMN matrix model [20j around a certain multi- fuzzy sphere solution is argued to be equivalent, 
using the Eguchi-Kawai reduction, to 4d N = 4 SYM. (For evidence supporting the validity 
of this formulation, see [21].) Given that the solution is BPS, it provides a regularization that 
preserve part of supersymmetry. Furthermore, given that one-dimensional system like the BMN 
model can be analyzed on computer, the previous formulation allows us to study 4d SYM using 
a Monte-Carlo simulation. 

A natural question to ask is to which kind of theories it is possible to apply the Eguchi-Kawai 
equivalence. The generalization to 4d N = 1 pure SYM is quite straightforward, as discussed in 
[22] . using the matrix model analyzed in [23]. However, introducing fundamental matter, as will 
be explained in § [3j turns out to be a difficult taslE In this paper, we show that quiver and Chern- 
Simons gauge theories can be regularized using the techniques of [16]. More specifically, in §[3] we 
construct the supersymmetric SU(N) x SU(M) gauge theory with bifundamental matter. Then, 
by sending the coupling constant of the latter gauge group to zero, we obtain a global flavor 
symmetry from this gauge symmetry, and as a consequence, this quiver gauge theory becomes 
supersymmetric QCD. In this construction both N c = N and Nf = MK (K is the number 
of bifundamental matters) must be infinite, but the ratio Nf /N c can be arbitrary. In § [J] we 
show that Chern-Simons theory, which is difficult to study on lattice, can also be formulated 
in terms of Eguchi-Kawai equivalence. For that purpose, we use a construction of the Chern- 

4 For other attempts that use the Eguchi-Kawai equivalence, see [171 1181 IT5] . 

5 As discussed in [24], large- N gauge theories with a quark in the two-index antisymmetric representation can 
be regarded as the counterpart of the usual QCD, in the sense that this representation reduces to the fundamental 
representation when N = 3. For these models, Eguchi-Kawai reduction can easily applied. See e.g. [15] , 
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Simons theory using a generalization of Taylor's T-duality prescription [26] which is discussed 
in [25]. Combining this results with the technique of [16] the Eguchi-Kawai formulation can be 
obtained straightforwardly. Furthermore combining the analysis for the quiver and Chern- 
Simons theories, we are able to construct the ABJM theory [TJ] from a matrix model. 

This paper is organized as follows. In § [2] we review the basic ideas of the Eguchi-Kawai 
equivalence. First, in § 12. 1[ we explain the quenched Eguchi-Kawai model \28\ 129] . Based on it, 
in § 12. 2\ we review the reduced model of SYM on S 3 [16] . In § [3] we generalize this technique 
to construct supersymmetric quiver gauge theories. In §[5] we formulate Chern-Simons theory in 
three dimensions along the line of [25] . Combining these results with the ones in § [3] we construct 
the Chern-Simons-matter theories which recently attracted much interest as the theory describing 
multiple M2-branes. 

2 The basics of the Eguchi-Kawai reduction 

In this section we review the Eguchi-Kawai equivalence [15j . In §2.11 we introduce the 
"quenched" version of the Eguchi-Kawai model [28^ I29| . which is relevant for our purpose. In 
§ 12.21 we use this technique to formulate large- X Yang-Mills on the three-sphere. 

2.1 Quenched Eguchi-Kawai model 

In the following we review the diagrammatic approach to the quenched Eguchi-Kawai model(QEK) 
|29| . The basic idea is that in the planar limit, Yang-Mills theory is equivalent to a matrix model 
around a suitable background. We will also consider QEK for compact space [16j [30] . In order to 
see clearly the difference between the compact and noncompact cases, we consider (analogously 
to [16]) the simplest example first, namely the correspondence between a zero-dimensional matrix 
model and a matrix quantum mechanics. 

As a simple example, we consider a matrix quantum mechanics with a mass term, 

Su = nJ dtTr Q(AX,) 2 - l[Xi, X,] 2 + ^Xf) , (1) 

where Xj (i = 1, 2, • • • , d) are N x N traceless Hermitian matrices. The covariant derivative Dt 
is given by 

D t X i = d t X i -i[A,X i }. (2) 
At large-X, this model can be reproduced starting from the zero-dimensional model 

Sbd = X ■ NTr H [y '^ ]2 " l [Xi > Xj]2 + IT*?) ' (3) 

where Y and Xj are N x N traceless Hermitian matrices. We embed the (regularized) translation 
generator into the matrix Y, 

Y b -9- = diag{p u - ■ ■ ,p N ), Pk = ±(k-?p\. (4) 

6 While this work was in progress we have been informed that the same idea had been studied by Ishiki et al 
[27] . We thank G. Ishiki for the discussion. 
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Figure 1: Two-loop planar and nonplanar diagrams with quartic interaction vertex. 



By expanding Y around this background, 

Y = Y b -9- + A, (5) 

the Feynman rules of the one-dimensional theory, as we will see in the following, are reproduced 
at large- N. 

The action can be rewritten as 



l i,3 



We add to it the gauge-fixing and Faddeev-Popov terms 

^. N Tr(±[Y b -9-,A} 2 -[Y b -9-,b][Y,c}) • (7) 

Then, the planar diagrams are the same as the ones in the Id theory up to a normalization factor. 
For example, a scalar two-loop planar diagram with quartic interaction (see FigJT]) is 

d(d-l)/l 2irN\ A (A/2vriV) (A/2vriV) 



- ) y — 

A / • ^Tli (Pi ~ Pk) 2 + m2 (Pj ~ Pk) 2 + m 2 

^ d(d-l) 2n f^ 2 dp f A / 2 dq 1 

4 A J_ A/2 2vr J_ A/2 2vr (p2 + m*)(q* + m*) " U 

The essence of this expression is that the adjoint action of the background matrix can be identified 
with the derivative and the matrix elements of the fluctuations can be identified with the Fourier 
modes in momentum space. The corresponding diagram in the Id theory is 

. vol . n 2 r ^ r ^ - 0) 

where Vol is volume of the spacetime. Hence by interpreting A and A/N to be UV and IR cutoffs, 
those diagrams agree up to the factor (y-) • Vol. The other planar diagrams also correspond up 
to the same factor. 
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The nonplanar diagrams do not have such a correspondence, but in an appropriate limit they 
are negligible. In the Id theory, by taking a planar limit they are suppressed by a factor N~ 2 . In 
the reduced model, they are suppressed if IR cutoff A/N goes to zero. To see this, let us calculate 
for example the two-loop nonplanar diagram in Fig{TJ It reads 

did- 1) A 

— is^s- (10) 

which is suppressed by a factor (A/TV) 2 compared with planar diagrams. 
Therefore, by taking the limit 

N^oo, A^oo, ( n ) 

the Id model on M. is reproduced from the Od model. 

If one wants to obtain the theory on a circle, it is necessary to fix the IR cutoff while suppressing 
nonplanar diagrams. This can be achieved by taking the background to be 

rb.g. ,. / ^ „ . A (, n l 



Y - g - = diagipi, ■ ■ ■ ,p ni ) ® 1„ 2 , p k = — (k - — , N = n\n<i, (12) 

n\\ I ) 

and taking m, n<i and A to be infinity while fixing the IR cutoff A/ni [16] 

It turns out that in this setup the background is not stable. So, to make the expansion 
meaningful, we have to "quench" the eigenvalues of Y, i.e. we have to fix the background by 
hand. This is the reason for the name "quenched" Eguchi-Kawai model. 

2.2 Eguchi-Kawai construction of Yang-Mills on S 3 

Next let us construct the Yang-Mills theory on the three-sphere by using the Eguchi-Kawai 
equivalence. The essence of QEK is to find a background whose adjoint action can be identified 
with the spacetime derivative. So, the strategy is to find a set of three matrices whose adjoint 
action can be identified with the derivative on S 3 . Such matrices were found in [3H 116] . In the 
following we will show the derivation in a heuristic way. 

2.2.1 YM on S 3 

In this section, we express the action of Yang-Mills theory on R x S 3 in a form convenient for 
our purpose [16]. The radius of the sphere is taken to be 2//i. The action of U(N) SYM is given 
by 

N 

J Js 3 

where is the 't Hooft coupling constant, g^ u (x) is the metric and g(x) is its determinant. The 
field strength is 



- J dt J d 3 xyf^o)Tr^F^, (13) 



= 8^ - d v A^ -i[A^ A v \. (14) 



7 Another reason to consider the ri2 — > oo limit is the following, if we take ni = 1 and quench the background, no 
zero-mode appears. This is not a problem when we take the noncompact limit, because the IR cutoff goes to zero. 
On the other hand, when we consider a compact space, the absence of the zero-mode destroys the Eguchi-Kawai 
equivalence. It turns out that by taking the 712 — > 00 limit, zero-modes are taken into account in an appropriate 
manner. 
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The Greek indices //, v refer to the Einstein frame and the Latin indices to the local Lorentz 
frame. 

The sphere part of this geometry has the group structure of SU(2). Given this group struc- 
ture, there exists a right-invariant 1-form dgg -1 and its dual Killing vectors £j, satisfying the 
commutation relation 

[Ci,Cj] = itijk£k- (15) 

Using the coordinates (9,ip,ip) defined by g = e -*¥" J 3/ 2 e - i6 '< :r 2/2 e -t^o-3/2^ ^ e vielbein E % can be 
expressed as 

E 1 = — (— sin ipdO + sin 9 cos ipdip) , (16) 
I 1 

E 2 = — (cos ipd9 + sin sin ipdip) , (17) 
A* 

£ 3 = - (dip + cos ddip) . (18) 



n 2 - 



In these coordinates the metric is given by 

dd 2 + sin 2 6 dtp 2 + (dip + cos 6 dtp) 2 . (19) 

The spin connection u a b c can be read off from the Maurer-Cartan equation, 

dE i - J jk E j A E k = 0, (20) 

/i 

^ijk = ~^ e ijk- (21) 



and the Killing vectors are given by 



d = —Efd M , (22) 
i 1 



Ci = 


—i ^— sin ipdg 


c 2 = 


—i ^COS (fdg — 


£3 = 


-id v . 



where 

3t 9 cos wag, H -oL 

sm # 

fl - r, . sin^ 

sm & 

(23) 

The Killing vectors represent a complete basis for the tangent space on S s . Furthermore given 
that the vielbeins are defined everywhere on S 3 , the indices i can be used as a label for the vector 
fields and l-formslf| 

By using the Killing vectors d, the action can be rewritten as [16] 



{f)^J dt J dShT 'fy D ' A '->' c - A)2 



This property is necessary in order to identify this index with the one in the matrix model [32| 
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+— (CiAj — CjAi) — —(CiAj — CjAi)[Ai, Aj] + ^[^4i 5 Aj] 
-^A\ + ifie^AiAjAk - i^AiiCjAk)) , 



(24) 



where Ai is defined in such a way that the 1-form of the gauge field on S 3 take the form A = AiE % 
and CZSI3 is the volume form of the unit three-sphere. 



2.2.2 Eguchi-Kawai reduction 

To construct matrices which represent derivatives on S 3 in a coordinate-independent way, it 
is useful to use the SU{2) group structure of S 3 . The Killing vectors (f23|) act on functions on 
S 3 ~ SU(2), whose irreducible decomposition ia^| 



C°°{SU{2))= Vj@---@Vj 

J=0,l/2,1,- y (2J+l)-times / 



(26) 



where Vj is the space that the spin J representation acts on. 

In order to realize this representation as the adjoint action of the background matrices, we 
first embed the SU (2) generators into N x N matrices. We then introduce the matrices Lj which 
satisfy the commutation relation of the SU(2) generators, 



[Li, Lj] — ieijkL 



k ■ 



(27) 



Since these matrices cannot be diagonalized simultaneously, we embed them in the following block 
diagonal form; 



\ 



L 



b's-1/2] 



L 



L 



S+1/2J 



V 



(28) 



where L^ s ^ is a (2j s + 1) x (2j s + 1) matrix which acts on the spin j s representation. The size of 
the matrix iV is 



N = J2(2j s + 1). 



(29) 



9 In general, for a compact Lie group G, a space of functions on G is decomposed as 

/ 



C°°(G) = 



V r • • • ffi Vr , 



(25) 



r \ d r —times 

where r runs over all the irreducible representations, V r is a representation space and d r is its dimension. 
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We introduce a regularization by restricting the representation space to a limited number of j s . 
Furthermore we take the integer s satisfying 

T T 

-2<*<2> (30) 
where T is an even integer. We introduce another integer P S> T and take j s to be 

P + s 



2 

The large N limit is taken in the following way 



(31) 



P^oo, T^oo, N^oo. (32) 

To see how this prescription works, let us consider the (j, j')-block, to which acts from left 
and L^'l acts from the right. A Basis for this block is symbolically written as 

\j,m)(j',m'\. (33) 

It can be decomposed into spin \j — j'\, ■ ■ ■ ,j+ j' representations. Let's count the number of 
representations of each spin. 

Spin : T + 1, because it appears only when j = j' > 0. 
Spin 1/2 : 2T, because it appears when j = j' ± 1/2. 

Spin 1 : (T + 1) + 2(T — 1) = 3T — 1, because it appears when j = f > 1 and j = j' ± 1. 
Spin JeZ: (T + 1) + J2i=i 2(T + 1 - 21) = (2 J + 1)T + 1 - 2J 2 . 
As long as T > J, we can approximate this expression as 

(number of spin J) ~ (2J + 1)T. (34) 

Therefore the representation space, or equivalently the variables appearing in the matrix model, 
can be regarded as a set of T copies of the space of functions on S 3 . As J increase the number of 
copies decreases. In this sense T plays a role of a momentum cutoff. 

In this way matrix elements can be identified with the functions on S 3 , or in other words the 
propagators in the Feynman diagram agrees. However it is not apparent if this identification is 
consistent with the multiplication of the fields. (This is necessary in order for the interaction 
vertices to agree.) In [16] it has been shown that 

| - (35) 

is a sufficient condition for the compatibility with the multiplication!^! 

10 Given that this condition force us to use very large matrix in a computer simulation, it would be nice if it 
could be relaxed. However this seems to be impossible because the eigenvalue distribution is not uniform without 
imposing T/P — > 0, while the eigenvalues should be distributed uniformly in order for the quenched Eguchi-Kawai 
equivalence to work. One possible solution is to make the density uniform by putting fuzzy spheres with the same 
radii on top of each other and tune the number of copies to be proportional to the spin. We thank G. Ishiki for 
stimulating discussion on this point. 
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By using these matrices we can relate a matrix model to a gauge theory on S 3 , given by the 
action (|24p . In order to do that, we consider the bosonic matrix quantum mechanics 

S = C 'idS ^ (\ {DtXi)2 + \ [Xh X]]2 + i ^ohX i X^X k - f^-Xf) , (36) 

where the constant C will be specified shortly. We then expand the action around a classical 
solution 



A t = 0, X t = -fiLt, (37) 
identify £j, A t and Ai with [Li, ■ ], A t and Xi + \iLi, 

Ci - [L h ■ }, Af d) -> A^ d \ A^Xi + pi*, (38) 
and replace the trace and the spatial integral by a trace, 



/' 



dn 3 Tr -> Tr. (39) 



The UV and IR momentum cutoff's are given by fiT and fi, respectively, and we will take the limit 
in such a way that 

\i -> 0, fiT -> oo. (40) 
In order to match the diagrams completely, the coupling constant should be taken as [16] 



16vr 2 , , 

A 4 d = Ai d 5-- 41 

In other words, we have to multiply the dimensionally reduced action by an overall factor 

C 1 ^. (42, 

This factor is analogous to the factor 2n/h in ([3]). Furthermore the four-dimensional 't Hooft 
coupling X^d should be scaled with the UV momentum cutoff \iT. 

Finally we would like to add a few remarks. First, the background is a classical solution and 
hence as long as it is stable we do not need to quench it. Second, when we take the large- N limit 
fixing the IR momentum cutoff fj,, in order to suppress the nonplanar diagrams it is necessary to 
change the background to 

- liLi ® l fc (43) 

and take k — > oo limit. Thirdly, this construction resembles the "twisted" Eguchi-Kawai model 
(TEK) [33] . In both cases the model is deformed by background flux terms so that noncommutative 
manifolds (fuzzy sphere for the former and fuzzy torus for the latter) become a classical solution, 
and the higher-dimensional theories are obtained as a fluctuation around these solutions. 
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2.3 Supersymmetry and stability of the background 

So far we have discussed only bosonic theories. Strictly speaking the Eguchi-Kawai equivalence 
does not work in bosonic models because the background is unstable. This problem arises quite 
generally, not only in the QEK [M], but also in the original reduction [28] and in the TEK 
[351 136| 157] . As usual, supersymmetry can remove such an instability [38] The advantage of the 
construction explained in § 12.21 is that supersymmetry can be preserved manifestly [16]. That is, 
the reduced model is a supersymmetric matrix model and the background (|28p is a BPS solution. 
Therefore we can expect the background to be stable at least at low temperature. 

For another approach to the Eguchi-Kawai reduction in supersymmetric Yang-Mills with uni- 
tary variables, see [18]. 

3 The Eguchi-Kawai model for quiver gauge theories 

In [22] 4d M = 1 SYM without flavor is formulated by using the Eguchi-Kawai construction 
of [16]. In order to consider QCD, we have to introduce flavors into this model. However it is 
difficult to describe fundamental matter along this line. The reason is the following. Because the 
derivative on the sphere is identified to the commutator with matrix i//Lj, the covariant derivative 
acting on the fundamental scalar ip can be written as 



The gauge field Ai acts only from the left, and Lj acting from the left can be identified with the 
background of corresponding field in matrix model, X{. However, the last term in the right hand 
side cannot be expressed as a matrix variable, since there are no field acting on ijj from the right. 
(In other words it does not appear from the dimensionally reduced model). To circumvent this 
problem, we consider bifundamental matter. Then the covariant derivative becomes 

Diip ~ i[/j,Li,ip] - iAiifj + iifjBi = i{[iLi - - if) ■ i(fJ.Li - Bi) = -iX^ + iipYi. (45) 

In this case, both Lj can be identified with the background of matrix variables and hence the 
technique of [16] can be applied. By taking the additional gauge coupling to be small, the gauge 
field Bi decouples and we restore the fundamental matter. 

3.1 Quiver matrix quantum mechanics and its quenched reduced model 

As the simplest example let us start by considering the bosonic quiver quantum mechanics 
with gauge group SU(N) x SU(M), where N and M are taken to be infinity by fixing the ratio 
M/N . For simplicity we take M = km and N = kn, where k, m, n are integers, and then, we take 
the k — > oo limit fixing m and n. We consider the following action 



Dii/j ~ i[/j,Li,ip] - iAiip = i(ixLi - Ai)ip - ip ■ i\iLi 



iXiip — if) ■ ijiLi. 



(44) 




(46) 



where (ft is N x M matrix and the covariant derivative acts on it as 



D t (p = d t (p - %A<\> + icpB. 



(47) 



ii 



An attempt to avoid the instability in the bosonic framework can be found in [39] , 
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Figure 2: A two-loop planar diagram of a bifundamental scalar <j). Solid and dotted lines represent 
SU(N) and SU(M) indices, respectively. 



Here A and B are gauge fields associated with SU{N) and SU(M), respectively. In this action, 
the field <j> is rescaled such that scaling parameter in k — > oo appears only in the overall factor, 
furthermore the parameters \i and g do not scale in this limit. 

This model is related to the reduced one via Eguchi-Kawai equivalence. The reduced model 
is given by 

^ • k Tr\-{X(f) - 4>Y){Y<jf - (f> ] X) + u 2 ^ + g{<j)(p ] ) 2 } . (48) 
If we expand this action around 

A / k\ 

X b - g - = diag(pi, ■ ■ ■ ,p k ) <g> l n , Y b - g - = diag(p 1 , ■ ■ ■ ,p k ) <g> l m , p r = — f r - - J , (49) 

then this model reproduce the results of the original one. 

As an example, consider the two-loop diagram shown in Fig. [2 In the matrix quantum 
mechanics, it gives 

VI * 2 f A/2 dp [V 2 dq 1 

Vol ■ k ■ on ml — / — — -= „. . n 50 

while in the reduced model it is 

2vr , 2 2 _ dp dq 1 



T' k - 9nm J^ /2 ^J^ /2 ^WT^WT7)- (51) 

Therefore we can see the correspondence as in § 12.11 up to the same factor (A/27r) • Vol. The 
generalization to other diagrams is straightforward. 

3.2 Bosonic quiver gauge theory in four dimensions 

Let us start by considering the bosonic quiver gauge theory with SU (N) x SU(M) gauge 
group. As in the previous subsection, we take the limit N, M — > oo with M = km, N = kn, 
k — > oo and m, n kept fixed. (As we will see, in terms of QCD with fundamental matter, this 
means N c , Nf — > oo with Nf/N c fixed.) Let us consider the action 



(52) 
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C-gauge — ^ 2 TfFpu ^2 ^ r ^ti" 

Cmatter = k (-Tr(D ^(D^rf + mjjTr^^j) , (54) 



where F^ u and G^ u are field strength of SU(N) and SU(M) gauge fields and B^, 

F„ v = d^A v - d v A„ - zL4 M , A,], G MV = d„B v - d u B^ - i[B„, B v ), (55) 

and 4>i(I = 1, • • • , if) are N x M complex matrices on which the covariant derivative acts as 

D^cpi = d^i - iA^dpi + tyiBy,. (56) 

We take the mass matrix mjj to be hermitian. 

In the Maurer-Cartan basis, the matter part of the Lagrangian reads 



^matter — kTv 



< (D t 4> I )(D t 4>i) J( - (ifiCifa - iAi4>i + i(t> I B i )(iiJ,Ci(j)i - iA^i + i4>iBrf + mjjTrcpjcf) 



The gauge part is the same as ([24)1 for each gauge group. 
By reducing the spatial dimensions to a point, we obtain 

C qX I qY i rtmatter /ro\ 

&id — b ld + i> ld + & ld , (be J 

S id = J J dtTr {^{DtXif + \\XuXjf + i^X^X* - f^Xfj , (59) 

SL = J dtTr Q(A*i) 2 + |[*i,^] 2 + ifieijkY i Y j Y k ~ ^ Y ^j > (60) 

gmatter = CkJ dtTr j (Atf/XA^j) 1 - - faYMXifa - faYtf + mijTrfatfj J ,(61) 

where Xj and 1^ are N x N and M x M scalar matrices which are obtained from A and -Bj. We 
take the background to be multiple of the background (|28p as in (I43p , 

X\- 9 - = -flLi ® l n , >f 9 ' = -M^i ® lm, (62) 

where the size of Li is k x k. Then, it is straightforward to see that the Eguchi-Kawai equivalence 
works around this vacuum. By taking g 2 B k to be zero, SU(M) reduces to (a part of) a global 
"flavor" symmetry. The number of flavors turns out to be Nf = MK, and the ratio Nf/N c can 
be arbitrary finite value. 

3.3 Supersymmetric quiver gauge theory in four dimensions 

We consider a quiver theory without superpotential. (incorporation of the superpotential is 
straightforward.) In this section we use the notation of Wess-Bagger's textbook [40] . 
The action is given by 

S = Sg aU g e + S ma tter = I dt I Xyg(x) gauge + ^-matter) , (63) 

J Js 3 
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Cgauge = Tr {--L (^F% + iXa^D^ - (^G% + iAs^A*) } , (64) 

kTrl-iD^D^rf - ifoiPD^i - iV24>\{\ A ^i - ^iX B ) 



-matter 



nq>j) 



(65) 



where ipj and <pi belong to (JV, M) representation as in the previous subsection. (Strictly speaking 
other multiplets are needed in order to cancel the gauge anomaly, but we omit them for notational 
simplicity. The modification is straightforward.) We notice that the last term is analogous to a 
mass term of adjoint scalars in 4d M = 4 on S 3 . 
The supersymmetry transformation is 

5^ A,j, = -iX A o^ + iea^XA, 

S^Bfj, = -iX B ane + iea u X B , 
S^Xb = G^e + ikgl^je, 
I = V2eipi, 



S^fa = i-fionDM + ^Zfa. (66) 

4V2 



Here the supersymmetry transformation parameter e satisfies 

D^e = ~^e. (67) 

The dimensionally reduced model is obtained by rewriting the action in the Maurer-Cartan 
basis and then by reducing the spatial dimensions. It is important that the parameters of the 
supersymmetry transformation depends only on t in this basis 

e(t) = e-^/ 4 e . (68) 

Consequently the dimensional reduction of spatial dimensions does not affect supersymmetry. 
By dimensionally reducing the spatial directions we obtain the matrix quantum mechanics 

L 7auge = jTr ^ {D t X t ) 2 + \{X h X,} 2 - ^X? + i^XiXjX^J 

C ( - 3 

+—Tr ( -iX x D t X x - X x o- l [Xi, X x ] - -fJ.X x X x 



9 A 



4' 



+ (X^Y), (69) 
LZ2er = CfcTr|(A0/)(A^) t -(X^ / -^)(X^ / -</» / y i ) t 
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-iV2(j)\(\xlpI - ^/Ay) + iV2(lpi\x - Xy^l 




The supersymmetry transformation reduces to 

5^X l = 

6^ A t = 

5^X X = 

S^Yi = 

sMBt = 

6^\ B = 



5^0! = 

S^i = 
where e is time-dependent 



i\xe + ieXx, 

-2(D t Xi) + e ijk [Xj, X k ] + 2%tiX^ a l e + ikg 2 A 4>i4>\^ 

iXyOie + iectiXy, 
iXye + ieAy, 

-2(D t Yi) + ^ k Y^Y k \ + 2^>y i ] a l e + ikgl^fae, 
V^a^X^i - faYi) + -^=e0/, 



4^/2 



e(t) 



e-^ 4 e . 



(70) 



(71) 



(72) 



It turns out that the background (|62|) preserve this supersymmetry. By expanding the ma- 
trix model around this background we recover the original 4d theory. Notice that we have to 
renormalize the bare 't Hooft couplings g\N and g B M appropriately in the continuum limit. 

Generalizations to more complicated quiver theories are straightforward. We emphasize that 
the equivalence works only when the field theory does not has gauge anomalies, because we 
assumed implicitly in the proof. 

The emergence of the chiral anomaly is a long standing problem in Eguchi-Kawai models. This 
problem exists in the present case too - the chiral symmetry seems kept in the reduced model. 
The chiral symmetry should be broken by some effect. Here, we do not pursue this direction, but 
just assume the presence of such an effect. One possible way to find it is to consider the chiral 
anomaly in the noncommutative space (concentric fuzzy spheres). 



4 The Eguchi-Kawai model for Chern-Simons gauge theories 

In this section, we consider the Eguchi-Kawai model of three dimensional Chern-Simons gauge 
theories. In a similar fashion to the YM cases, we can obtain matrix models which have a fuzzy 
sphere as a classical solution. Quivers can be introduced into this construction as in the previous 
section. We consider the ABJM model j!4j as an example. 
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4.1 Bosonic case 

Let us start with the bosonic Chern-Simons with U(N) gauge group. The action is 

Scs = i-^j d 3 x^{x)e^Tr (-A^d u A p - jA^A^j , (73) 

where k is an integer. On the three-sphere of radius 2/fi, by using the Maurer-Cartan basis, this 
action can be written as [25] 

Scs = i ~ J d 3 nTr e^ k {-(C t A 3 )A k + {C 3 A t )A k ) - jAAjA^ + pA? . 

(74) 

Dimensionally reducing it, we obtain 

k rrl ( 2 Z mic^j- -wr T 7- """2 



S mm = iC ■ —Tr I — — € iwp X ll X v Xp + fiXf j . (75) 

This theory has a classical solution 

Xi = -fiLi, [Li,Lj] = ie ijk L k . (76) 



It is easy to check that the 3d action ([?4"|) is obtained from ([75]) by taking the background ([28]) 
and using the mapping rule given in § 12.21 

4.2 ABJM theory 

The previous construction can be easily promoted to supersymmetric Chern-Simons theory. 
As a concrete example, let us formulate the ABJM model, which gives a description of the multiple 
M2-branes theory [TJ] (see also [H]). In the following we take the gauge group to be U(N) x U (M). 
When this model is put on the three sphere, an additional mass term must be added in order to 
keep super symmetry. The model is given by 

ik ( 2i 2i 

c = —^PTr\-A li d v A p --ArA v A p + B li d v B p + -B ll B v B p 

+^-Tr (D^D^a + i^D^ a ) + kV(cP, VO, (77) 

where 

~Tr {{ct> a 4> a f + $ a <f> a ) 3 ) - ^Tr («^^> 7 ^) 
la j. lax 13\ , V 



+-Tr [M^pf^) + ^<f> a <t>«- (78) 
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Here and are gauge fields of U(N) and U(M) groups, respectively. Bifundamental matter 
fields 4> a ,ip a (a = 1, ■•• ,4) are in the (N,M) representations and the covariant derivative 
acts as 



D 



V 



iAa<t> a + i4> a Bn, 



(79) 



and similarly for ip a . We have rescaled these matter fields in such a way that the Chern-Simons 
level k appears only in the overall factor. 
The supersymmetry transformation is 



ri^o-^^a + r) a f 3 a IJL (j) a 'if) 13 

2 
3 



P, 



'-] 



r? a + - 



(80) 



where the parameter satisfies rj a/3 = —rf a , (Tj a p)* = \if L ^ri^ = if 1 ^ and V^r/"^ = — ^a^"q a/3 . 

Rewriting the action using the Maurer-Cartan basis and taking the zero-dimensional reduction, 
we obtain the matrix model 



>5,, 



Ck 



i f 1i 2? 
—Tr l-jeWXiXjX* + fiXf + je^Y^Y, - ^Y? 



-^-Tr(Y^ a - ^ a Xi)(X^ a - (paY 1 ) 



(81) 



Expanding this action around the background (|28p . we can reproduce the original action. As 
before, the transformation parameter in the 3d theory is a constant in the Maurer-Cartan basis 
and the reduced model is supersymmetric. The supersymmetry transformation is 









- (V^ 


s (od) Y . = 


- (V^ 




a\-iX, 



^iP + 2 -e a ^^^r, Sp - ^W 7 . 

(82) 



Notice that this background preserves all supersymmetries. 

In the actions the Chern-Simons level k appears as an overall factor. The original ABJM is 
reproduced from the reduced model in the planar limit with both k/N and k/M kept fixed. 
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5 Conclusion and Discussions 

In this paper we have applied a recently proposed large-iV reduction technique [16] to su- 
persymmetric quiver and Chern-Simons theories. As concrete examples we have constructed the 
SU(N) x SU(M) supersymmetric quiver gauge theory with bifundamental matter fields and the 
ABJM model of multiple M2 branes. Furthermore, by taking one of the gauge couplings to be 
small in the supersymmetric quiver gauge theory we obtain SU(N) supersymmetric QCD with 
flavor. In this construction both N c and Nf are infinite but the ratio Nf/N c can take any value. 
Therefore this construction provides us with a valuable tool to study the dynamics of supersym- 
metric QCD, e.g. supersymmetry breaking, Seiberg duality conjecture [32], etc. 

The reduced model of the ABJM theory will be useful to study the AdS/CFT correspondence 
numerically. Of particular interest is the strong 't Hooft coupling region that is expected to 
describe type IIA string on AdS 4 x CP 3 . This region can be studied by using the Eguchi-Kawai 
equivalence. It turns out that the parameter region where k is smaller than O(N) is also important 
to obtain insights into M-theory. Thermo dynamical properties are also interesting. For these 
reasons it is still valuable to study a lattice formulation which is valid at finite- N and can be put 
at finite temperature. For references in this direction, see e.g. [43] . 

We expect that these models have the sign problem, unlikely to the reduced model for 4d 
J\f = 1 pure SYM [22] . This problem possibly make it difficult to study these models numerically. 
At finite temperature, however, the sign problem might be mild, similarly to the case of the 
maximally supersymmetric matrix quantum mechanics [10]. It is important to study by direct 
simulation how severe the sign problem is. 

Although inherently restricted to the planar limit, the Eguchi-Kawai equivalence can be a 
powerful tool to explore the dynamics of supersymmetric gauge theories. Numerical studies on 
these models will be reported in future communications. 
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